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rfy Abstract. General Relativity (GR) describes gravitation well at the energy scales 

which we have so far been able to achieve or detect. However, we do not know whether 
GR is behind the physics governing stronger gravitational field regimes, such as near 
neutron stars or massive black-holes (MBHs). Gravitational-wave (GW) astronomy is 
a promising tool to test and validate GR and/or potential alternative theories of grav- 
ity. The information that a GW waveform carries not only will allow us to map the 
strong gravitational field of its source, but also determine the theory of gravity ruling 
5_i \ its dynamics. In this work, we explore the extent to which we could distinguish be- 

tween GR and other theories of gravity through the detection of low-frequency GWs 
from extreme-mass-ratio inspirals (EMRIs) and, in particular, we focus on dynamical 
Chern-Simons modified gravity (DCSMG). To that end, we develop a framework that 
enables us, for the first time, to perform a parameter estimation analysis for EMRIs in 
■s^j- " DCSMG. Our model is described by a 15-dimensional parameter space, that includes 

(f-) . the Chern-Simons (CS) parameter which characterises the deviation between the two 

■ theories, and our analysis is based on Fisher information matrix techniques together 

£N| ' with a (maximum-mismatch) criterion to assess the validity of our results. In our anal- 

ysis, we study a 5-dimensional parameter space, finding that a GW detector like the 
Laser Interferometer Space Antenna (LISA) or eLISA (evolved LISA) should be able 
| to discriminate between GR and DCSMG with fractional errors below 5%, and hence 

place bounds four orders of magnitude better than current Solar System bounds. 
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1. Introduction 



EMRIs are low-frequency GW systems detectable by GW observatories like (e)LISA. 
They are made up of a MBH (M. ~ 10 4 - 10 7 M Q ) and a stellar compact-object (SCO) 
(m + « 1 - 50M Q ). Once the SCO is captured by the MBH, the system loses energy 
and angular momentum and the SCO starts inspiralling driven by the emission of GWs. 
During this process, the geometry of the MBH spacetime where the inspiral take place 
is mapped in the shape of the GWs emitted. Therefore, by looking at the phase and 
timing of an EMRI waveform, we will be able to determine the physical parameters of 
the system and the structure of the MBH potential. For this reason, waveform models 
used for GW detection rely on the structure of the MBH spacetime and, hence, on the 
theory of gravity assumed. 

GR has been well tested in the weak field, but the behaviour of gravity in these 
weak gravitational regimes could be effectively described by GR, whereas in strong 
gravitational fields the true theory might deviate from GR. In that case, gravity may 
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be described as an effective theory whose action contains higher-order curvature terms. 
Then, since EMRIs emit deep inside the strong MBH potential, through the detection 
of their GWs we can unveil the physics behind these systems and reveal the theory of 

gravity that governs them. 

In this contribution, based on ref. dCanizares et al.1 12012 ). we focus on study- 
ing the detectability of DCSMG, a modification of GR constructed by the addition 
of a Chern-Simons (CS) gravitational term (also known as the Pontryagin invariant) 
to the action. Interest in this theory was initiated with the work of Jackiw and Pi 
ref. dJackiw & Pill2003h where gravitational parity violation was investigated. In ad- 
dition, DCSMG presents distinctive features that could be observed by GW detection, 
such as higher-order curvature terms which beco me important in s trong gravitational 
fields. This work follows the conventions of ref. dMisner et al.ll 1973b . 



2. Parameter Estimation in DCSMG 

In order to study EMRIs in DCSMG, we need to model the orbit of the SCO and the 
GWs generated by the system in this theory, which means that we have to know the 
metric of the spacetime where the SCO moves. Besides from the black hole mass 
M. and spin a, black hole metrics in DCSMG are parameterized by a single universal 
constant, £ := a 2 /(fiK N ), where k n is the gravitational constant, and a and ft are universal 
coupling constants that control the strength of the CS modifications. Then, £ is the only 
parameter responsible for the deviations between GR and DCSMG. The non-vanishing 
components of the DCSMG metric are ref. (lYunes & Pretoriusll2009l) : 

/ 2M.r\ p 2 2 X 2 



5 £ a Mi 12M. 27M 2 \ 2M,ar 
r- 1 + + 



8M?M. r 4 \ 7r 10r 2 / p 2 



sin 2 8 , 



with p 2 = r 2 + a 2 cos 2 8, A = r 2 f + a 2 , f - 1 - 2MJr, and I = (r 2 + a 2 ) 2 - a 2 A sin 2 8. 

Once we have the metric of the M BH spacetime in DCSMG, we obtain the geodesic 
equations ref. (|Yunes & Pretorius1l2009h : 

i = i K + L z 5gf{r), 4> = 4> K -E6gf(r), (2) 

r 2 = r 2 K + 2ELJ6gf(r) 8 2 = 8 2 K , 

where the dots denote differentiation with respect to proper time and E and L z are 
respectively the energy and the angular momentum per unit SCO mass. The quantities 
(i K , r K , 8 K , K ) represent the geodesic equations in GR, that is in Kerr, and the CS 
deviations are governed by the term 6g^ s <x g ■ a ref. ( Sopuerta & Yunesll2009b . 



Since the orbit of the SCO is not a geodesic trajectory, but an inspiral that evolves 
adiabatically, the orbital parameters associated with the geodesic should be evolved to 
model the trajectory. To that end, we have developed an sche me to evolve the EMRI or- 
bit in the framework of DCSMG ref. (|Canizares et alJl2012h by incorporating radiation 
reaction (RR) effects, which not only drive the inspiral, but also break the degeneracy 
between an orbit in GR and another in DCSMG. In other words, the RR pushes the SCO 
across the MBH geometry, which differs in both theories, and hen ce, the GW gener- 
ated in GR is dephased with respect to the one generated in DCSMG ( Sopuert a" & Yunes! 
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2009h . Once the orbital evolution of the SCO is obtained, the waveform associated with 
its tra jectory is co mputed through a multipolar expansion up to the mass quadrupole 
ref. dThornejll98dh . 




Figure 1. Symmetric normalized overlap [Eq. ©], which shows the projection of 
a GW in DCSMG with the corresponding one in GR. 

We describe EMRI systems with a 15 -dimensional parameter space including the 
CS parameter A = [M„ a, p, e Q , p Q , # inc q, f = g ■ a,9 s ,(f> s ,9 K ,(/> K ,D L ,i// ,x ,<f> } (see 
Table [J for definitions), where the subscript refers to values taken at the initial time 
of the inspiral and ip and x are angle variables (AV) introduced to describe the ra- 
dial and polar motion respectively. However, the present work is restricted to a 5- 
dimensional subset of the parameters that show the strongest correlations with £,. The 
reduced parameter space is then given by A = {M„a,e ,£,D L }. Here, we assume 
a/M, = e Q = 0.25, and £/M^ = 5 • 10~ 2 and the rest of the parameters have fixed values 
given by: p = mJM, = 10 M Q /M„ p = 11, 6» mc0 - 0.569, 6 S = 1.57, S - 1.57, 
K = 0.329 , D L = 1 Gpc, and (f> K = ^ = Xq = (f> = 0.78 . 



Parameter 


Description 


Parameter 


Description 


M. 


MBH mass [M e ]. 


s s 


EMRI polar angle. 


a = \S.\IM. 


MBH Spin [Af.]. 




EMRI polar angle. 


li = mJM. 


EMRI mass ratio. 




MBH spin polar angle. 


e o 


Eccentricity. 




MBH spin azimuthal angle. 


Pq 


Semilatus rectum. 


D h 


Distance to the EMRI [Gpc]. 


^inc.O 


Orbit inclination . 




AV. for the radial motion. 




Xo 


AV. for the polar motion. 


c 


?-a [Mil 


^0 


BL azimuthal angle. 



Table 1. Summary of the parameters that characterize an EMRI system in DC- 
SMG. Units are indicated in square brackets for parameters with dimension and Di 
is measured from the SSB to the EMRI. 



In our analysis, we consider two different EMRI system, labelled A and B, char- 
acterised by the mass of their MBH. In particular, for system A we consider M. = 
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5 ■ 1O 5 M and M m = 10 6 M Q for system B. The parameter estimation errors expected 
for LISA for these systems are given in Table|2] 

The signal-to-noise (SNR) for eLISA is around two times smaller than the one of 
LISA, then, to compare the error estimates of both detectors we can normalize to a fixed 
SNR. The results obtained for system A show that the parameter estimation accuracy 
does not change appreciably when the noise curve of LISA is changed for the one of 
eLISA and so, all previous results can be considered to apply to eLISA as well, with 



the corresponding SNR corrections ref. (Can izares et al.1 12012). 



System A 


System B 


AlogM. ~ 5 ■ 10 _a 
Aa ~ 5 • 1(T 6 M. 

Ae () ~ 3 • 10- 7 
Alog f ~4 - 1(T 2 
Alog(D L /^)~2-10- 2 


AlogM. ~ 6- 1(T 4 
Aa ~ 3 • 1(T 6 M. 

Ae ~ 10" 7 

Alog f ~2 - 1(T 2 

Alog(D L / / /)~2-10- 2 



Table 2. Parameter estimation errors for LISA observations of the systems. 

In order to study how different a waveform in DCSMG, h cs , is from the corre- 
sponding one in GR, h GR , we can compute their (normalized) projection, i.e. their 
symmetric normalised overlap function: 

u 1 ( h GRl h cs) . 

O [hGR, h cs j ^ - = , (3) 

Vl h GRl h GRA h Csl h Csj 

where both waveforms correspond to EMRI systems defined with the same parameters. 

In Figure 1, we show how the projection of h cs onto h GR changes as the values 
of the MBH spin a/M, and the CS parameter £/M^ are modified. These results are 
for system A, assuming a 0.5yr EMRI observation. We see that, for higher values of 
a/M, and f/M? the overlap Eq. © decreases, since the difference in the evolution of 
the SCO in G R and CS increases and , consequently, the deviations of h cs from h GR are 
enhanced ref. dCanizares et alj|2012h. 

In our work ref. dCanizares et al.ll2012h . we have performed numerical simulations 
to explore how the results change under different assumptions, namely the effects of (i) 
including RR, (ii) changing the value of the spin and (ii) changing the value of the CS 
parameter. The parameter es timation anal ysis is performed using the Fisher information 
matrix (FM) formalism ref. ( Fishej[T935h . and in order to assess the validity of the FM 



results, we evaluate the maximum-mismatch criterion (MMC) ref. (IVallisneril 12008). 
Regarding (i), we evolved system A with and without RR and for different evolution 
times: r gvol = 0.1 , 0.3 , 0.5 , 1 yr, which ensures that enough high SNR is accumulated 
to perform parameter estimation. We conclude that the inclusion of the RR improves 
the SNR of the signals and also the parameter estimates, in particular those of the spin, 
a, and of the CS parameter, This is partially due to the increase of the overall SNR 
due to RR, but even after rescaling to a fixed reference SNR, we see an improvement in 
the parameter measurement accuracies when RR is included. As one could expect due 
to the adiabatic nature of the RR, the improvement with the inclusion of RR is more 
significant for longer evolution times. On the other hand, to explore (ii), that is the 
effect of the MBH spin a on the error estimates, we consider that systems A and B have 
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the following values for the MBH spin: a/M. - 0.1, a/M. - 0.25, and a/M. - 0.5. 
We found that the smaller the spin, a/M., the better the parameter estimate for the 
CS parameter, £. In particular A£ ~ 2.8 • 10~ 2 for system A and A£ ~ 1.4 • 10 
for system B. The reason for this is that the CS modifications affect a single MBH 
metric component, g^ [Eq. ©], which contains the CS parameter f/M., multiplied 
by the spin parameter a. The unperturbed Kerr metric component is proportional to a, 
and so the relative change in this metric coefficient due to the addition of the DCSMG 
correction is proportional to Since we keep £ = a £ fixed as we vary a, the value 
of £, increases as a decreases and so the CS correction to the MBH metric is larger 
relative to the leading order Kerr metric term. Finally, we address (iii) by varying the 
value of the CS parameter £/M 4 of system A. In particular we consider: £ = 0.05M 4 , 
£ = 0.1M 4 and £ = 0.2M?, finding that the value of the CS parameter £/M 4 only affects 
significantly the error estimate of the CS parameter itself, whereas modifying a/M. has 
a major effect on the error estimates of all the parameters employed in our study, and in 
particular on £. 

With the framework that we have developed, we can go an step further and assess 
whether LISA could bound the value of the CS parameter by estimating the error on 
the measurements of A£, i.e. by setting a bound like £ < A£, assuming that GR is 
the true theory of gravity. As different EMRI systems will provide different constraints, 
and £ is a universal quantity, we just need to look for the EMRI system that provides 
the best constraint. In this way, we find that (e)LISA could place a bound < 
1.4 • 1 4 km, which is almost four orders of magnitude better than th e bound obtained 
in ref. dAli-Haimoud & Chenll201lh using Solar System data (see ref. dYagi et alj|2012l) 
for a recent work where the CS parameter is estimated using BH binaries). 

In summary, we have analysed how well a space-based GW detector like (e)LISA 
can discriminate between an EMRI system in GR and one occurring in a modified 
gravity theory like DCSMG. To that end, we have constructed a waveform template 
model in DCSMG and used the Fisher matrix formalism, together with the MMC, to 
estimate errors in parameter measurements, finding that a GW detector like (e)LISA 
should be able to discriminate between the two theories with a fractional error below 
5%. Furthermore, assuming that GR is the true theory of gravity, we have been able to 
constrain the CS parameter £, with bounds four orders of magnitude better than using 
Solar System data. 
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